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ABSTRACT 


The population structure, particularly growth, age, mortality and somatic production of the 
olivid snail Olivancillaria deshayesiana were investigated. Annual formation of internal shell 
growth marks was confirmed by the record of stable oxygen isotopes in the shell, which reflects 
seasonal patterns of water temperature. A von Bertalanffy growth model fitted to 81 size-at- 
age data pairs, indicating that O. deshayesiana may attain 31 mm SL in about 10 years. The 
estimated total mortality Z and natural mortality M were 0.651 y-1 and 0.361 y-1, respectively. 
Fishing mortality F was 0.290 y-1, and the exploitation rate E was 0.445, indicating that this 
population was not overexploited at the time of the study. However, this situation may well 
change in the future, since the important prawn and shrimp fisheries (in intensity and scale) 
in the Mar del Plata area (38°20’S, 57°37’W) may indirectly affect the exploitation status of 


the studied population. 
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INTRODUCTION 


The determination of the individual age within 
a population is important for a comprehensive 
understanding of population dynamics and 
is required to establish growth and mortality 
patterns. Several methods have been used to 
study growth in gastropods (e.g., Shigemiya & 
Kato, 2001; Richardson et al., 2005a, b; Bar- 
roso et al., 2005), including the family Olividae 
(Stohler, 1969; Tursch et al., 1995; Rocha- 
Barreira, 2003; Caetano et al., 2003). In recent 
years, the technique by Epstein et al. (1953) 
based on the analysis of the isotope compo- 
sition of calcium carbonate shells has been 
increasingly used by many authors (Giménez 
et al., 2004; Cledon et al., 2005; Bigatti et al., 
2007; Arrighetti et al., 2010). 

The family Olividae is composed of marine 
gastropods of moderate size, rarely exceed- 
ing 10 cm in total shell length. They inhabit 
Sandy coastal shores in tropical and temperate 
regions of the world (Smith, 1998). There are 
eight Olivancillaria species in the southwest- 
ern Atlantic Ocean (Teso & Pastorino, 2011), 
some of which are the target of subsistence 


fisheries (Scelzo et al., 2002; Narvarte, 2006). 
This family, together with some volutids (i.e., 
Adelomelon, Zidona) and nassarids (i.e., Bucci- 
nanops) represent the most important endemic 
gastropods of sandy bottoms in the Argentine 
malacological province. 

Olivancillaria deshayesiana is a small gas- 
tropod (up to 35 mm) distributed from Rio de 
Janeiro (22°59’S, 43°11’W), Brazil, to Neco- 
chea (38°33’S, 58°44’W), Argentina, at 0—60 
m depth (Teso & Penchaszadeh, 2009; Teso 
& Pastorino, 2011). Individuals are semi-buried 
in the substratum, and their shell is partially 
covered by the mantle. The surface of the shell 
exhibits a distinct pattern of growth bands that 
have been used for ageing purposes, as in O. 
auricularia (Rocha-Barreira, 2003). In addition, 
Caetano et al. (2003) estimated growth and 
mortality of O. vesica from length-frequency 
data with computer-based methods. 

Little is known about the reproduction and 
ecology of O. deshayesiana. Borzone (1995) 
described its embryonic development and egg 
capsules, and Teso & Penchaszadeh (2009) 
reported the imposex phenomenon in a popula- 
tion from the Mar del Plata area. This species 
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is taken as bycatch in the fisheries for prawn 
Artemesia longinaris and shrimp Pleoticus 
muelleri, rather than being the direct target of 
artisanal fishermen (Scelzo et al., 2002). The 
objective of this study was to determine growth, 
mortality and the individual production of O. 
deshayesiana in the Mar del Plata area and to 
evaluate its significance for the local fishery. 


MATERIALS AND METHODS 


Specimens of Olivancillaria deshayesiana 
were collected monthly by bottom trawl (3 m 
mouth opening, 15 mm of mesh size) off Mar 
del Plata City (88°20’S, 57°37’W) between 5 
and 12 m depth between October 2005 and 
September 2007. In this area, water tem- 
perature ranges between 8.4 and 20.8°C and 
salinity is rather constant at 33.8%o (http://www. 
lafe.uba.ar/tele/Antares_argentina/; Guerrero & 
Piola, 1997). Animals were taken to the labora- 
tory, where the shell was separated from the 
soft body. Shell length (SL) was measured with 
a vernier caliper to the nearest 1mm and spiral 
growth trajectory length (SG) with a thread to 
the nearest 1 mm (Fig. 1). Shell-free wet mass 
(SFWM) was weighed to the nearest 0.1 g. 
The relationship between SL and SFWM was 
calculated. 


Age and Growth 


Only the pattern of internal growth marks was 
considered for the analysis due to the presence 
of columellar material forming a callus on the 
spire. The shells of ten individuals were cut 
along the whorls following the spiral growth 
trajectory from the apex to the posterior end of 
the aperture (Fig. 1), as described by Bigatti et 
al. (2007) and Arrighetti et al. (2010), among oth- 
ers. The cut surface was ground with sandpaper, 
using grits of P400, P1200, P2400 and P4000 
grade. The polished shell cuts were examined by 
stereomicroscope. Visible internal growth marks 
were identified and the corresponding growth 
trajectory length SG; from the apex to each mark 
i was measured as described by Cledon et al. 
(2005) and Gurney et al. (2005), among others. 
A stable oxygen isotope ratio (8180) was used 
to confirm the annual periodicity in the forma- 
tion of growth marks. 5180 in biogenic calcium 
carbonate structures is inversely proportional 
to temperature during carbonate deposition, as 
demonstrated empirically by Epstein & Lowen- 
stam (1953) and Epstein et al. (1953). 

Two representative specimens of 31 mm and 
29 mm SL were used for stable isotope analysis. 
The carbonate samples (+ 50 ug each) were 
collected from the central shell layer along the 
SG by means of a dental drill (bit size 0.5 mm) 


FIG. 1. Olivancillaria deshayesiana. FIG. 1A, B: Showing the spiral growth trajectory (SG), along which 
the samples for oxygen and carbon isotopes were taken, dash line in A indicates location of section in 
B; FIG. 1C: Detail of the internal growth marks (arrows). Scale bar = 1 cm. 
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at a sample-to-sample distance of about 1.5 
mm. Stable oxygen isotope composition was 
determined in the stable isotope laboratory of 
the Alfred Wegener Institute with a Finnigan 
MAT251 mass spectrometer coupled to an 
automatic carbonate preparation device. The 
precision of measurements was better than + 
0.08% for 618O based on repeated analysis of 
a laboratory-working standard over a one-year 
period. Data are related to the Pee Dee belem- 
nite (PDB) standard through repeated analyses 
of National Bureau of Standard (NBS) isotopic 
reference material 19 by: 

om [(Reample / Retandard) z 1] * 1000 

where Rsampie and Retandara are the isotopic re- 
lationship of the sample and of the standard, re- 
spectively. 618O (d18O,) values were converted 
into water temperature by the paleotemperature 
equation of Epstein et al. (1953), as modified by 
Craig (1965) and Anderson & Arthur (1983): 

T(°C) = 16.00 — 4.140 * (d180, — d18O,,) + 

0.130 * (d18O, — d180,,)2 

Sea water 8180 (d'8O,,)) was estimated at 
-0.8%o from the gridded data (version 1.1) of 
LeGrande & Schmidt (2006); see also the 
corresponding NASA website (http://data.giss. 
nasa.gov/o18data). Visual shell marks coincid- 
ing with 618O peaks were interpreted as winter 
growth marks, and the distance between two 
subsequent winter growth marks as the annual 
shell increment. Thus, each winter growth mark 
i and its corresponding shell growth trajectory 
length SG; were interpreted as one size-at-age 
data pair. SG; was translated into shell length 
SL; by means of the linear relationship between 
the two parameters. The iterative non-linear 
Newton algorithm was used to fit a growth 
model to the size-at-age data. The specialized 
von Bertalanffy growth model was found to fit 
the data best: 

SL,= SL. * (1 — e-K*(t-to)) 

where SL. is asymptotic shell length, K is 
the body growth coefficient, t is age and tọ is 
a theoretical age at which length equals zero. 
Overall growth performance was estimated by 
the index: 

P = log(K * SL..3) 

as described by Moreau et al. (1986) & Pauly 
et al. (1994). 


Production 


Individual somatic production was calculated 
by the mass-specific growth rate method ac- 


cording to Brey (2001) from the size frequency 
distribution, the size-body function and the 
growth model parameters: 

P;= BM; * G; 

where BM, is the mean individual body mass 
at size j and G is the annual mass specific 
growth rate at size j as calculated by: 

GORE SL/SM 

Smaller size classes are likely to be under- 
sampled due to the selectivity of the com- 
mercial bottom trawl, and thus population 
production (P) was not calculated. 


Mortality, Yield per Recruit and Exploitation 
Rate 


It was presumed that mortality in O. 
deshayesiana is adequately described by the 
single negative exponential mortality model and 
the corresponding size-converted catch curve 
(Pauly, 1984a, b) was used to estimate total 
mortality rate Z (y-1), 

N;/At; = No * gai 

where N; is number in size class j, Afjis time 
required to grow through size class j, and t, is 
age at midsize of size class j. Natural mortality 
rate M was estimated by means of an empirical 
model that relates production-to-biomass ratio 
(P/B) ratio to maximum age fina, (y), maxi- 
mum body mass BMmax (kJ) and mean water 
temperature T (Kelvin) in unexploited benthic 
populations (Brey, 1999, 2001). If growth can 
be described by the von Bertalanffy model 
and mortality by the single negative exponen- 
tial model, mortality rate M and P/B ratio are 
equivalent (Allen, 1971), hence: 

log(M) = 1.672 + 0.993 * log(1/fmax) — 0.035 * 
log(BM ax) — 300.447/T 

N = 901, £ = 0.880 

Maximum age and body mass were inferred 
from the largest animal through the growth 
model and the size-mass relationship (applying 
the conversion factor 3.818 kJ/gWM, average 
for Gastropoda, Brey, 2001). 

Fishing mortality F and exploitation rate E 
were computed by: 

F =Z-M, andE=F/Z 

A selection curve giving the probability of 
capture PC; per size class j was computed by 
a procedure outlined by Pauly (1984a, b). Yield 
per recruitment was computed by the Beverton 
& Holt (1964) model, from which a predicted 
value for optimum exploitation rate was also 
obtained. 
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FIG. 2. Size-frequency distribution of Olivancillaria deshayesiana (N 
= 420) from Mar del Plata area in 2005-2007. Superimposed dotted 
line indicates individual somatic production Pi (g shell-free wet mass 
year-1) per size class. 
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FIG. 3. 6180 profile along a transect from apex to shell edge from two 
individuals of Olivancillaria deshayesiana. FIG. 3A: 31 mm SL, 10 y; FIG. 
3B: 29 mm SL, 7 y. Arrows indicate position of internal growth marks. 
Water temperature (right-hand scale) computed from carbonate 18O by the 
paleotemperature equation of Epstein et al. (1953); see text for details. 
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FIG. 4. Von Bertalanffy growth curve for Olivancil- 
laria deshayesiana SL,= 38 * (1 — e-0.147*(t-3.392})); r 
2 = 0.913; N= 81 size-at-age data pairs obtained 
from the 10 specimens analyzed. 


RESULTS 


A total of 470 snails were collected, which 
ranged from 16 to 35 mm in total shell length 
(Fig. 2). Shell-free wet body mass ranged from 
0.4 to 10.8 g and was related to SL by: 

log(SFWM) = -4.76+ 3.60 *log(SL); (g - mm; 
2 = 0.77; N = 470) 

SL was linearly related to SG, the distance 
from the apex along the growth trajectory to the 
aperture of the shell, by: 

SL = -42.87 + 3.54 * SG; (mm — mm; £ = 0.93; 
N = 30) 


Age and Growth 


The 6180 profiles showed distinctly oscillat- 
ing pattern in both shells (Fig. 3). The average 
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FIG. 5. Size-converted catch curve based on the 
size-frequency distribution (Fig. 2) and growth func- 
tion (Fig. 4) of Olivancillaria deshayesiana. Dots: 
data included in regression; circles: data excluded 
from regression. In (N/dt) = -0.651 ti + 8.193; 7 
= 0.93; 95% confidence range of slope = + 0.164. 


range of -1.1%o to + 0.8% 180 corresponded to 
a temperature range of about 9.8°C to 17.4°C. 
Almost all 8180 maxima (corresponding to tem- 
perature minima) coincided with a shell growth 
mark, thereby indicating that one growth mark 
is formed each winter. 

The von Bertalanffy growth model: 

SE = 38 (=e aE. N= 81, r2 = 
0.913 

fitted best the 81 size-at-age data pairs ob- 
tained from the 10 specimens analyzed (Fig. 4). 
Overall growth performance of O. deshayesi- 
ana was P= 3.91. 


Production 


Individual somatic production P; increased 
with size j up to a maximum of 0.116 g SFWM 
year? at 25 mm SL and decreased towards 
larger sizes (Fig. 2). 


Mortality, Yield per Recruit and Exploitation 
Rate 


Total mortality rate Z of O. deshayesiana 
amounted to 0.651 y-1 (95% confidence range + 
0.164, Fig. 5). Natural mortality rate M was esti- 
mated at 0.361 y-1, with tnax = 20.7 y (at 35 mm 
SL), BMmax = 41.23 kJ (10.8g SFWM) and T = 
286 K. Accordingly, fishing mortality F amounted 
to 0.290 y-1 and exploitation rate E to 0.445. Fig- 
ure 6A shows the selection curve for the currently 
applied sampling technique. Specimens less 
than 17 mm were not caught at all; probability of 
capture attained 50% at 22.8 mm SL and 100% 
2 27 mm SL. The optimum exploitation rate was 
estimated at 0.770 (Fig. 6B). 


DISCUSSION 


Stable Oxygen Isotope Ratios and Environ- 
mental Temperature 


In the shell of O. deshayesiana, stable oxygen 
isotope ratios reflect seasonal patterns of water 
temperature; each maximum 8180 value corre- 
sponds to a winter minimum temperature. Thus, 
the distinct shell mark at each 180 maximum 
indicates a slowing shell growth during winter, 
as observed for other gastropods from the Mar 
del Plata region (Cledon et al., 2005; Gimenez 
et al., 2004). The average annual temperature 
inferred from the 8180 profile in the shell ranged 
between 9.8°C and 17.4°C, in agreement with 
the reported seasonal instrumental tempera- 
tures of 8.4—20.8°C in the same area (http:// 
www.iafe.uba.ar/tele/Antares_argentina/). 
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FIG. 6. Olivancillaria deshayesiana. FIG. 6A: Plot 
of probability of capture versus length in Olivan- 
cillaria deshayesiana; FIG. 6B: Yield per recruit- 
ment in Olivancillaria deshayesiana based on the 
Beverton and Holt (1964) model (calculations by 
means of the FISAT software reference). 


Age and Growth 


From the isotopic data and the internal growth 
marks, it is inferred that O. deshayesiana 
may reach 31 mm SL in about 10 years. This 
maximum age is similar to those reported for 
Olivella biplicata (8—12 years) from California, 
USA (Stohler, 1969) and Oliva oliva (10 years) 
from Papua, New Guinea (Tursch et al., 1995), 
while it is distinctly above those estimated for 
O. auricularia (6 years) from Rio Grande, Brazil 
(Rocha-Barreira, 2003) and O. vesica (4—5 
years) from Rio de Janeiro, Brazil (Caetano 


et al., 2003). In general, cold-water marine 
invertebrates found at higher latitudes grow 
more slowly and live longer than closely related 
species occurring in warmer waters (Dehnel, 
1955; Ray, 1960; Fonseca et al., 2000). This 
may account for the shorter life span of O. 
vesica compared with species of the same 
genus living at higher latitudes. 

The specialized von Bertalanffy model de- 
scribes the shell growth of O. deshayesiana 
properly, because it is a gastropod in which 
shell growth has no inflexion point (e.g., Frank, 
1965; Caetano et al., 2003; llano et al., 2004; 
Richardson et al., 2005a; Chatzinikolaou & 
Richardson, 2008; Eversole et al., 2008; Ar- 
righetti et al., 2011). This is in contrast to other 
gastropods with a sigmoid pattern of shell 
growth (e.g., Giménez et al., 2004; Cledon et 
al., 2005; Bigatti et al., 2007). Therefore, all 
standard fishery analysis based on the von 
Bertalanffy model, such as virtual population 
analysis, applies to O. deshayesiana. 

The overall growth performance of O. 
deshayesiana (P = 3.91) is below the range 
found for other large commercially important 
gastropod species from temperate regions 
(reviewed in Arrighetti et al., 2011). 


Production, Mortality and Exploitation Rate 


The size-frequency distribution was based 
on data from catches with nets of 15 mm mesh 
size, and therefore sampling efficiency was 
lower than 50% for individuals smaller than 
23 mm SL (Fig. 6A). This is clearly seen in the 
size-converted catch curve, where individuals 
younger than 5 years (27 mm SL) are under- 
represented. 

Therefore, estimates of production and 
productivity based on this distribution would 
distinctly underestimate true population values. 
The estimated fishing mortality of 0.290 y-1 and 
the current exploitation rate of 0.445 indicate 
that this population was not overexploited at the 
time of the study. However, this situation may 
well change in the future, since the important 
prawn and shrimp fisheries (in intensity and 
scale) in the Mar del Plata area (Scelzo et al., 
2002) may indirectly affect the exploitation sta- 
tus of the studied population. Furthermore, O. 
deshayesiana is a food source for commercial 
fishes such as Urophycis brasiliensis (Brazil- 
ian codling), Sympterigia bonapartei (ray), 
Callorynchus callorynchus (elephant fish) and 
Micropogonias furnieri (white croaker) (Pen- 
chaszadeh et al., 2006). It is possible that the 
depletion of the population may cause a shift 
to a new stable trophic situation. 
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